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THE  SACLANTCEN  OCEANOGRAPHIC  DATA  BASE 
Vol .  II:  ACCESS,  INTERROGATION  AND  DISPLAY 


by 

Richard  F.J.  Winterburn 


ABSTRACT 


An  oceanographic  data  base  established  at  SACLANTCEN  on  an  'in-house' 
UNIVAC  1106  computer  system  is  described.  Volume  I  discusses  the  design 
criteria  used  in  setting  up  the  data  base,  lists  its  structure  and  content, 
and  explains  how  acquired  data,  either  from  outside  institutions  or  from 
SACLANTCEN  experiments  are  re- formatted  and  entered.  Volume  II  describes 
how  data  are  accessed,  interrogated,  and  displayed,  including  the  plotting 
of  charts  with  coastlines  and  of  contoured  data. 

4 

INTRODUCTION 


An  oceanographic  data  base  has  been  established  using  an  in-house  UNIVAC 
1106  computer  system.  Volume  I  <1>  describes  the  background  to  the 
establishment  of  thj  data  base,  discusses  the  design  criteria  used  in 
setting  it  up,  lists  its  structure  and  content,  and  explains  how  acquired 
data  are  reformatted  and  entered. 

The  present  document,  Volume  II,  goes  on  from  that  point  by  describing  the 
methods  and  facilities  developed  to  access,  interrogate,  and  display  the 
contents  of  the  data  base.  The  suite  of  programs  described  has  been 
written  to  produce  "standard"  outputs  of  profile  display,  spatial  and 
temporal  distribution  charts,  and  computations  of  mean  profile  and 
gradient.  A  simple  interface,  also  described,  may  be  included  in  any  user's 
FORTRAN  program  to  make  the  contents  of  the  data  base  available  for  other 
purposes. 

The  executive  programs  are  designed  primarily  to  be  used  in  demand  mode 
from  interactive  graphics  terminals.  They  call  on  a  library  of  routines  to 
select,  display  and  analyse  data  from  the  base,  using  as  far  as  possible 
fully  conversational  programming  techniques.  This  almost  entirely  divorces 
the  user  from  the  data-base  structure  and  also  from  the  application 
programs;  thus  if  one  of  the  outputs  described  in  this  report  is  suitable 
for  a  users'  needs  he  may  produce  it  with  no  knowledge  of  any  programming 
language  or  data  management  system.  The  only  prerequisite  is  an  apprecia¬ 
tion  of  the  oceanographic  content  and  a  basic  knowledge  of  English. 

The  function  and  use  of  each  of  those  routines  is  described  in  Ch.  5, 
together  with  examples  of  the  output.  A  summary  of  the  software  require¬ 
ments  is  given  in  Ch.  1,  followed  by  a  description  of  the  more  important 
supporting  facilities.  The  successful  development  of  the  executive  routines 
depends  on  efficient  data-field  definition  and  retrieval,  and  well- 
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supported  graphic,  cartographic,  and  interactive- terminal  software;  these 
are  described  in  Chs.  2,  3  and  4,  respectively. 


1  SOFTWARE  REQUIREMENTS 

As  the  data  base  has  been  created  for  a  number  of  well-defined  applications 
its  retrieval  capabilities  should  closely  reflect  the  output  requirements 
of  those  applications.  These  may  be  summarized  as  follows: 

□  To  make  meso-scale  and  large-scale  climatic  studies  of  oceanographic 
features  in  SACLANTCEN' s  areas  of  interest. 

□  To  provide  an  integrated  environmental  data  system  for  the  use  of 
those  SACLANTCEN  projects  that  require  pre-cruise  appraisal  and 
post-cruise  analysis  of  oceanographic  conditions. 

□  To  allow  research  into  the  methodology  for  an  historical  oceano¬ 
graphic  data  base  with  a  view  to  improving  its  statistical 
description  for  use  as  input  to  oceanographic  forecasting  and 
acoustij  prediction  models. 

a  To  provide  historical  ocean-surface  data  to  be  used  for  inter¬ 
comparison  with  remotely-sensed  data  and  real-time  surface  data. 

These  applications  therefore  require  the  retrieval  of  data  both  in  space 
and  time  and  the  description  of  data  individually,  summarily,  and 
collectively. 


2  DATA- RETRIEVAL  FIELD  DEFINITION 

2. 1  General 

Every  retrieval  of  data  from  the  data  base  requires  a  clear  field  defini¬ 
tion;  this  is  done  in  conversational  mode  by  a  routine  called  CAMPO  (see 
Sect  2.5). 

As  was  described  in  detail  in  Vol.  I  <1>,  the  data  are  stored  by  means  of  a 
system  of  attribute  keys  and  within  every  data  record  a  dictionary  area 
contains  additional  unique  information  about  that  particular  profile,  each 
item  of  which  may  be  regarded  as  a  key.  These  keys  can  be  thought  of  as 
primary  and  secondary  in  the  sense  that 

(a)  Primary  keys  are  those  that  provide  direct  access  to  a  data  record 
by  virtue  of  its  file  name/element  name, 

Instrument  -  Marsden  Square  (MSQ)  -  One-degree  Square  (DSQ)  -  Month 

(b)  Secondary  keys  are  those  that  are  stored  within  the  data  record  but 
may  be  used  for  oata  selection, 

Position  -  Date  -  Time  -  Country  -  Ship  -  Cruise 

These,  however,  may  only  be  interrogated  sequentially  within  one  series  of 
primary  key  combinations. 
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Data  retrieval  is  therefore  facilitated  by  the  use  of  these  keys,  and  may 
be  conveniently  divided  into  three  distinct  data-field  definition  criteria 
i.e.  Spatial,  Temporal,  and  Source. 


2,2  Spatial  Definition 

Spatial  definition  is  possible  by  using  one  of  the  four  methods  displayed 
graphically  in  Fig.  1. 


2.2.1  Marsder,  area 

This  allows  retrieval  of  data  from  a  particular  Marsden  square  or  one  of 
its  1°  sub-squares  (Fig.  2).  It  is  the  quickest,  most  direct  of  the  four 
methods  as  no  computation  is  required  to  create  the  data-file  keys,  whereas 
the  other  methods  require  additional  input  data  and  computation  with  which 
to  create  the  keys. 


2.2.2  Rectangular  area 

This  allows  retrieval  within  a  given  geographical  area,  whose  limits  are 
defined  as  shown  in  Fig.  3a. 

These  limits  may  be  input  directly  as  signed  numerical  values,  or,  they  may 
be  input  interactively  using  the  cursor  of  a  graphic-display  terminal.  The 
latter  method  is  defined  as  shown  in  Fig.  3b  when  input  on  a  chart 
displayed  on  the  terminal  screen.  From  these  cursor  values  the  limits  are 
computed  as  values  of  latitude  and  longitude,  .hich  then  allows  the 
computation  of  the  MSQ/DSQ  coverage. 


2.2.3  Circular  Area 

This  allows  retrieval  within  a  given  circular  area  defined  by  the  position 
at  the  centre  and  a  range  in  nautical  miles  (Fig.  4).  The  geographic 
limits  of  the  north,  south,  east  and  west  extremities  of  this  circle  are 
computed,  followed  by  the  relevant  MSQ/DSQ  coverage. 


2.2.4  Ground  Track 

This  allows  retrieval  within  an  area  defined  as  a  swath  along  a  computed 
ground  track  between  two  positions  as  defined  in  Fig.  5  .  The  swath  width 
is  the  total  width  in  n.mi.  From  those  input  data  the  geographic  extrem¬ 
ities  of  the  catchment  area  are  computed,  followed  by  the  MSQ/DSQ  coverage. 

During  the  data- retrieval  phase,  the  data  are  retrieved  using  the  MSQ/DSQ 
coverage  but  only  those  data  falling  within  the  relevant  area,  be  it 
rectangular,  circular  or  a  track  swath,  are  extracted  from  the  data  base. 
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1»  SQUARE  IDENTIFICATION 
NORTHERN  HEMISPHERE 


10°  WEST  0°  EAST  10° 


10°  0°  10° 


FIG.  2  MARSDEN  SQUARE  IDENTIFICATION 


FIG.  3(a)  RECTANGULAR  AREA  DEFINITION-DIRECT 


FIG.  3(b)  RECTANGULAR  AREA  DEFINITION-INTERACTIVE 
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FIG.  4  CIRCULAR  AREA  DEFINITION 
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Long 
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2.3  Temporal  Definition 

The  month  o  observation  is  a  direct  key  by  which  data  may  be  retrieved 
without  examining  the  record  contents.  Within  the  dictionary  area,  Year, 
Day  and  time  of  day  are  recorded  and  may  be  used  for  selection  purposes; 
the  standard  field  definition  software,  however,  offers  only  Month  and  Year 
as  a  selection  key. 


2.4  Source  Definition 

Within  this  group  are  four  parameters  that  may  be  defined  for  data 
retri eval 

(a)  Instrument  (b)  Ship  (c)  Country  (d)  Cruise  (SACLANTCEN  only). 

Of  which  (a)  is  a  direct  key,  (b)  and  (c)  are  dictionary  items  as  described 
in  <1>,  «nd  (d)  is  an  in-house  identifier  used  only  for  data  recorded  by 
SACLANTCEN. 


2.5.  Conversational  Selection  Routine  (CAMPO) 

This  routine  is  called  by  each  executive  program  to  define  the  data 
selection  within  the  fields  described  previously.  It  is  totally  conversa¬ 
tional  by  requesting  each  data  field  to  be  defined,  using  a  philosophy  of: 

HOW  MANY?  and  WHICH? 

a  response  of  999  to  the  HOW  MANY?  question  signifies  that  all  available 
data  within  the  field  are  required,  in  which  case  the  second  question  is 
skipped. 

All  responses  are  solicited  by  the  SOE  (>)  character. 

It  should  be  noted  that  in  all  data  descriptions  within  this  memorandum, 
north  latitude  and  east  longitude  should  be  input  as  positive,  whereas 
south  latitude  and  west  longitude  should  be  input  as  negative  values. 

The  demands  made  by  this  routine  are: 

HOW  MANY  INSTRUMENTS  DO  YOU  WISH  TO  SEARCH 
ANSWER  1  TO  5  OR  999  FOR  ALL 

> 

PLEASE  LIST  THE  INSTRUMENT  CODES,  WHERE 

1  FOR  MBT  *  NON- SACLANTCEN 

2  FOR  XBT  *  NON- SACLANTCEN 

3  FOR  NC  *  NON- SACLANTCEN 

4  FOR  SACLANTCEN  XBT  DATA 

5  FOR  SACLANTCEN  STD  DATA 
> 

YOUR  DATA  AREA  MAY  BE  DEFINED  SPATIALLY 
BY  USING  ONE  OF  THE  FOLLOWING  METHODS: 

1  BY  DEFINING  A  MARDSEN  SQUARE/DEGREE  SQUARE  LIST 

2  BY  DEFINING  NORTH,  SOUTH,  EAST  AND  WEST  GEOGRAPHICAL  LIMITS 

3  BY  INTERACTIVE  MAP  SELECTION  (DEMAND  MODE  ONLY) 

4  BY  DISTANCE  FROM  A  POSITION 

(I.E.  WITHIN  A  CIRCLE  OF  RADIUS  DISTANCE) 

5  BY  DEFINING  IN  A  TRACK  AND  SWATH 
> 
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If  the  response  is  1: 

A  MAP  OF  THE  ATLANTIC  AND  MEDITERRANEAN  COASTLINES 
WILL  NOW  BE  PLOTTED.  WHEN  THE  CURSOR  IS  ENABLED, 
INPUT  THE  SOUTH-WEST  AND  NORTH-EAST  CORNERS 


This  will  be  followed  by  a  report  on  the  area  selected  and  the  MSQ/DSQ  list 
will  be  computed. 

If  the  response  is  2: 

PLEASE  INPUT  YOUR  AREA  LIMITS  IN 

DEGREES  AND  MINUTES  OF  LATITUDE  AND 

LONGITUDE  IN  THE  ORDER  TOP,  BOTTOM,  RIGHT  AND  LEFT 


The  MSQ/DSQ  list  will  now  be  computed. 

If  the  response  is  3: 

HOW  MANY  MARSDEN  SQUARES  DO  YOU  WISH  TO  SEARCH? 
ANSWER  EITHER  1  TO  23 
OR  999  FOR  ALL 

> 

PLEASE  LIST  THE  MSQ  NUMBERS 

> 

FOR  MARSDEN  SQUARE  [repeated  for  each  MSQ  requested] 
HOW  MANY  DEGREE  SQUARES  DO  YOU  WISH  TO  SEARCH? 

ANSWER  EITHER  1  TO  99 
OR  999  FOR  ALL 

> 

PLEASE  LIST  DEG.  SQU.  NUMBERS 


If  the  response  is  4: 

PLEASE  INPUT  THE  POSITION  IN  TERMS  OF 

LATITUDE  AND  LONGITUDE  (DEGS  AND  MINS  S  AND  W  NEGATIVE) 

FOLLOWED  BY  THE  DISTANCE  IN  N.MI. 

If  the  response  is  5: 

PLEASE  INPUT  THE  TRACK  AND  POSITIONS 
IN  TERMS  OF  LATITUDE  AND  LONGITUDE 
DEGREES  (I)  AND  MINUTES  (F) 

SOUTH  AND  WEST  NEGATIVE 
FOLLOWED  BY  THE  TOTAL  SWATH  WIDTH 
IN  N.MI. 

> 

MANY  YEARS  DO  YOU  WISH  TO  SEARCH? 

ANSWER  EITHER  1  TO  70 
OR  999  FOR  ALL  AVAILABLE 

> 

PLEASE  LIST  THE  YEARS  (THE  LAST  TWO  DIGITS  E.G.  71  FOR  1971) 

> 
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HOW  MANY  MONTHS  DO  YOU  WISH  TO  SEARCH? 

ANSWER  EITHER  1  TO  11 
OR  999  FOR  ALL  AVAILABLE 

> 

PLEASE  LIST  THE  MONTHS  (NUMERIC  E.G.  1  FOR  JANUARY) 

> 

IN  HOW  MANY  COUNTRIES  DO  YOU  WISH  TO  SEARCH? 

ANSWER  EITHER  1  TO  20 
OR  999  FOR  ALL  AVAILABLE 

> 

PLEASE  LIST  THE  COUNTRY  CODES  (E.G.  35  FOR  FRANCE) 

> 

FOR  HOW  MANY  SHIPS  DO  YOU  WISH  TO  SEARCH? 

ANSWER  EITHER  1  TO  20 
OR  999  FOR  ALL  AVAILABLE 

> 

PLEASE  LIST  THE  SHIP  CODES  (ALPHA  E.G.  CA  FOR  CALYPSO) 

> 


3  CARTOGRAPHIC  DATA 
3.1  General 

In  the  selection  and  display  of  oceanographic  data  it  is  obvious  that 
cartography  plays  an  important  role  in  clarifying  the  spatial  charact'  — 
i sties  of  the  data. 

For  data  selection,  the  interactive  method  (Sect.  2.2.2)  requires  a 
graphic- termi nal  display  of  coastline.  For  data  display  in  mapping  and 
horizontal  contouring,  again  a  plotted  coastline  is  a  valuable  enhancement 
to  the  output.  For  these  reasons  a  number  of  data  sets  have  been  created  at 
SACLANTCEN  by  digitizing  the  coastline  in  particular  areas  of  interest; 
these  are  described  in  Sect.  3.2  Horizontal  distribution  charts  and  contour 
plots  are  produced  on  the  standard  Mercator  projection,  which  is  the  most 
widely  used  projection  in  navigation.  This  is  briefly  explained  in 
Sect.  3.3,  together  with  details  of  the  geodetic  spheroid  options 
avai Table. 


3.2  Coastline  Digitization 

For  a  small-scale,  low-resolution  representation,  the  entire  seaboard  of 
the  N.  Atlantic  Ocean  has  been  digitized  at  a  scale  of  1:39  000  000,  as 
shown  in  Fig.  6;  this  is  used  as  a  general-purpose  coastline.  As  shown  in 
Fig.  7,  when  a  small  part  of  this  is  selected  for  enlargment  the  low 
resolution  of  the  digitization  becomes  apparent.  For  this  reason  two 
further  data  sets,  see  Figs.  9  and  10.,  have  been  created  for  areas  of 
particular  interest.  These  have  been  digitized  at  a  much  larger  scale  and 
the  increased  detail  is  apparent  by  comparing  Fig.  7  with  Fig.  8  which  is 
the  same  area  but  plotted  with  data  digitized  at  1:9  000  000. 
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These  data  sets  are  selected  and  accessed  by  routines  COAST  and  COSTA 
respectively,  the  former  soliciting  the  input  parameters  required  by  the 
latter,  as  follows: 

PLEASE  INDICATE  FROM  WHICH  COASTLINE  DATA  SET 
YOU  WISH  TO  PLOT: 

-1  NO  COASTLINE 

1  MEDITERRANEAN  (1:12M) 

2  S.W.  APPROACHES  TO  ENGLISH  CHANNEL  (1:9M) 

3  ATLANTIC  OCEAN  (1:39M) 

> 


If  the  response  is  positive  then 

PLEASE  SELECT  THE  LINE  WIDTH  TO  BE  USED: 

1  THIN  LINE 

2  THICK  LINE 

3  DASHED  LINE 

> 


3.3  Chart  Projection  and  Geodetic  Data 

All  charts  produced  by  this  suite  of  programs  use  the  standard  Mercator 
projection,  which  is  satisfactory  within  the  latitudes  of  our  areas  of 
research.  At  higher  latitudes,  however,  a  conical  projection  is  more 
appropriate. 

Pearson  <2>  defines  the  full  numerical  solution  and  gives  an  excellent 
overview  of  the  advantages  and  disadvantages  of  these  and  other 
projections.  The  mapping  programs  however  make  use  of  a  "graphical" 
solution  by  computing  the  projection  dimensions  through  Meridional  parts, 
which  is  the  normal  navigation  computation  as  used  by  Hydrographic  Offices 
<3>. 

The  Meridional  parts  are  computed  by  using  routine  MERPAR,  which  allows  a 
choice  of  spheroidal  constants  to  define  the  reference  ellipsoid  between 
the  two  most  widely  used  Geoids,  viz.  the  International,  or  Hayfords,  1910 
and  the  World  Geodetic  System  (WGS)  1972,  the  primary  constants  for  which 
are  defined  by  the  semi -major  axis  a  and  the  flattening  f, 

where  f  =  - — -,  in  which  3=  semi -mi nor  axis, 
a 

The  relevant  constants  being: 

International 
a  -  6378388.0  metres 
f  =  1/297 

WGS 

a  =6378135.0  metres 
f  =  1/298.26 
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4  SUPPORTING  SOFTWARE 

4. 1  Data  Retrieval 

An  in-house-developed  data  management  system  <4> ,  known  as  the  Mini -Filing 
System  (MFS),  is  used  for  all  data  retrieval  and  allows  direct  access  to 
individual  data  records.  The  routines  of  this  processor,  which  are  used  by 
the  FINDBT  routine  to  extract  individual  profiles  are  not  documented  here. 
The  FINDBY  routine  is  therefore  the  interface  between  all  the  executive 
programs  and  the  data  base  by  using  the  MFS  routines 

OPRDF 

REDES 

REDATA 

CLOS 

SETMD 


4. 2  Graphics 

An  in-house-developed  graphics  system  <5>  is  used  to  create  device¬ 
independent  files  of  plot  instructions,  i.e.  the  files  may  be  subsequently 
output  to  either  Calcomp  plotter  or  Tektronix  screen.  The  "cal comp- like" 
calls  of  this  system  are  used  in  many  routines  of  the  SMODS  package  but  are 
not  documented  here.  This  package  is  available  world-wide  as  it  is  now 
included  in  the  UNIVAC  users'  library. 


4.3  Contouring 

A  complete  package  of  contouring  and  three-dimensional  plotting  routines 
<6>  has  been  implemented  at  SACLANTCEN.  This  extremely  powerful  package 
(CCP)  allows  gridding,  smoothing  and  contouring  of  irregularly  spaced  data. 
It  has  been  adapted  and  used  by  the  CONTOUR  and  SECTION  executive  programs 
through  routines 

INSIDE 

ZGRID 

SMOOTH 

GETLEV 

CONSEG 

DATAPT 


5  DISPLAY  AND  ANALYSIS  SOFTWARE 


5. 1  General 

The  application  software  consists  of  six  executive  routines  that  call  on  a 
library  of  44  sub-routines  and  also  on  the  previously  mentioned  UNI*TEKX, 
MFS  and  CCP  packages.  The  main  programs,  written  in  FORTRAN  V  are  known  as 

DISPLAY 

MAP 

MEAN 

CONTOUR 

SECTION 

SIGMAT 
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and  are  described  in  App.  A.  The  function  of  each  of  the  subroutines  is 
described  in  App.  B  and  the  links  between  executive  programs  and 
sub-programs  are  illustrated  graphically  with  each  main  program  in  Figs. 
A.l,  A. 3,  A. 5,  A. 7,  A. 9,  and  A. 11. 

Each  executive  routine  is  described  in  terms  of  function,  execution, 
input-data  format,  and  output-data  examples.  It  will  become  apparent  that, 
as  far  as  possible,  conversational  mode  programming  is  employed,  almost  all 
the  input  being  solicited  by  plain-language  (English)  requests.  In  addition 
all  responses  are  checked  by  the  I/O  routines  for  "reasonableness"  i.e. 
invalid  responses  are  rejected  and  the  particular  request  is  repeated  but 
borderline  cases  can  pass  as  good  input  data. 


CONCLUSIONS 

The  use  of  a  suite  of  interactive  programs  to  access,  interrogate  and 
display  data  of  the  SACLANTCEN  Oceanographic  Data  Base  has  been  described. 

These  programs  produce  "standardised"  outputs  of  profile  display, 
horizontal  and  vertical  distribution,  mean  profile  computation  and  density 
displays. 

A  simple  interface  has  been  developed  that  allows  any  FORTRAN  program  to 
access  the  data  within  the  data  base. 

As  the  data  base  is  expanded  and  its  application  increased,  further  reports 
will  be  produced  to  update  this  users'  description. 
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BATCH 


APPENDIX  A 
EXECUTIVE  PROGRAMS 


A. 1  Display  Program 

Function;  To  display  either  graphically  and/or  numerically,  single  or 
multiple  profile-  of  temperature,  salinity  or  sound  speed  as  a  function  of 
depth;  or  to  compute  for  any  profile,  values  of  potential  temperature, 
<A.1>,  sigma-t  <A.2>,  density  <A.2>,  potential  density,  sound  speed  <A.5>, 
and  Brunt-Viisala  frequency  <A.4>. 

Execution: 

aXQT  SMODS*SYSTEM. DISPLAY 
Subroutines:  (see  Fig.  A.l) 


DISPLAY 


t . .  r  - 

SETMD  BIGPIC 


SITETl) 


1 

HEAD 

1 - 1 - 1 — 

CAMPO  NBGA  EXPLAN 

- 1 - 

INCHEK 

1 - 1 - 

TESTA  MNDBT 

- 1 - 1 - 

LIST  CALC 

- T" 

UKLIST 

r 

UPLOT 

.J 

1 

BIST 

| 

STANO  NTS 

1 

TESTA 

L 

1 

TESTA 

1 

TESTA 

i  1 

MERPAR 

r- 

•  I 

JANGLE 

ALLSQ 

LIMITS  UN1*TEKX 

VBU 

Vic  WILSV 

ROVO 

PRES 

SIG  POT 

TPLOT 


NBGN 


OUTPUT 

I 

NBGN 


UNI*TEKX  GENMAR 


I  I  I  |  I  I 

PLOBOR  ALNER  COSTA  MERPAR  LONLEN  UNI*TEKX 


UN1*TEKX  MERPAR  J 

MERPAR 

MERPAR 


L 


CCP 


FIG.  A.l  ROUTINES  OF  EXECUTIVE  PROGRAM  -  DISPLAY 


Input  Data 

Field  Definition  via  subroutine  CAMPO  (see  Sect.  2.5  of  main  text)  followed 
by: 


WHICH  TYPE  OF  OUTPUT  DO  YOU  REQUIRE? 

TYPE  1  LISTING  ONLY 

TYPE  2  CALCOMP/TEKTRONIX  PLOT  ONLY 

TYPE  3  LISTING  AND  CALC/TEKT  PLOT 

TYPE  4  COPY  TO  MT 

TYPE  5  COPY  ON  DISK  VIA  MFS 


SACLANTCEN  SM-lOi 


If  a  listing  is  required 

WHICH  STYLE  OF  LISTING  DO  YOU  REQUIRE? 

STYLE  1  HEADING  ONLY 

STYLE  2  HEADING  PLUS  OBSERVED  LEVELS 

STYLE  3  HEADING  PLUS  STANDARD  LEVELS 

STYLE  4  HEADING.  OBS.  AND  STAND.  LEVELS 

STYLE  5  COMPUTATIONAL  LISTING  AT  OBSERVED  LEVELS 

STYLE  6  COMPUTATIONAL  LISTING  AT  STANDARD  LEVELS 


NOTE:  If  run  in  batch  mode  or  on  the  Tektronix  4014  the  print  output  has  a 
character  plot  of  each  trace. 

If  a  plot  is  required: 

WHICH  PARAMETER  DO  YOU  WISH  TO  PLOT? 

TYPE  1  FOR  DEPTH  vs  TEMPERATURE 
TYPE  2  FOR  DEPTH  vs  SALINITY 
TYPE  3  FOR  DEPTH  vs  SOUND  SPEED 

> 

WHICH  TYPE  OF  PLOT  DO  YOU  REQUIRE? 

TYPE  1  SINGLE  PROFILE  OF  EACH  STATION 
TYPE  2  MULTIPLE  PROFILES  (ENVELOPE) 

> 

PLEASE  TYPE  MIN  AND  MAX  PARAMETER  AXIS  VALUES 


PLEASE  TYPE  MIN  AND  MAX  DEPTH  AXIS  VALUES 


PLEASE  TYPE  SCALE  FACTOR  OF  PLOT 
WHERE  1.0  GIVES  A  PLOT  5  INS  BY  8  INS. 

> 


If  a  copy  on  M.T,  is  required 

PLEASE  ENSURE  A  TAPE  HANDLER 
IS  AVAILABLE  TO  BE 
ASSIGNED  TO  YOUR  RUN 
IF  ALL  OK  TYPE  "GO" 


If  a  copy  on  disk  via  MFS  is  required: 
PLEASE  TYPE 

QUALIFIER  FILENAME  (12  CHARS.  EACH) 
IN  WHICH  RECORDS  SHOULD  BE  WRITTEN 


The  input-data  solicitation  is  terminated  by: 

THIS  CONCLUDES  YOUR  OUTPUT  OPTIONS 
IF  ALL  OK  TYPE  999 
ELSE  TYPE-  999 
> 
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If  the  print  option  has  been  specified  and  the  run  is  being  executed  in 
demand  mode  the  following  request  is  made 

00  YOU  WANT  YOUR  OUTPUT  SENT  TO  THE  LINEPRINTER? 

If  the  response  is  positive,  the  data  is  buffered  into  a  temporary  mass- 
storage  file  using  the  name  21,  which  must  be  output  before  the 

a  FIN  control  statement. 

If  the  response  is  negative,  the  output  is  directed  to  the  terminal  screen 
and,  if  needed,  a  hard  copy  may  be  taken  with  the  4631  unit. 

Output  Data:  (see  Fig.  A. 2) 
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FIG.  A. 2 (a)  DISPLAY  -  DESCRIPTIVE  AREA  PRINTOUT 
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FIG.  A. 2(b)  DISPLAY  -  PRINTOUT  OF  ALL  DEPTH  LEVELS 
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explanation  of  column  headings 
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137.500 

.006 

-.001 

16.371 

-97.593 

-.001 

-91.942 

175.000 

.007 

.001 

15.867 

84. 782 

.001 

89.430 

2P5.000 

.006 

.000 

16,124 

236.55? 

.000 

953.781 

275.000 

.006 

.000 

16.130 

232.058 

.000 

1068.370 

350.000 

.006 

.001 

16.118 

183.496 

.000 

268.742 

450.000 

.006 

.001 

16.134 

182,658 

.000 

269.184 

580.000 

.006 

-.000 

16.252 

-235,389 

-.001 

-167.681 

650.000 

.007 

.001 

16,027 

101.724 

.001 

1)1.173 

750.000 

.006 

.000 

16. 22? 

296.506 

-.or.,-. 

-489.460 

850.000 

.007 

.001 

16,098 

113.250 

.001 

1  ?5,  845 

950.000 

.007 

.001 

16.097 

107.419 

.on 

113. 722 

1050.000 

.006 

.001 

16.156 

122.135 

.001 

129.281 

1150.000 

.006 

.001 

16.217 

1 45,544 

.001 

153.632 

FIG.  A. 2(c)  DISPLAY  -  COMPUTATIONAL  LISTING  AT  STANDARD  LEVELS 
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DEPTH  t  fl  J 


DEPTH  ! 
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A. 2  MAP  Program 

Function:  To  plot  a  Mercator  chart  including  intermediate  lines  of  latitude 
and  longitude,  intermediate  border  tick  marks,  a  coastline  and,  if, 
required  select  and  plot  data  from  the  data  base. 

Execution: 

a  XQT  SMODS*SYSTEM. MAP 
Subroutines:  (see  Fig.  A. 3) 


PLOBOR  ALMER  COSTA  MERPAR  LONLEN  UNl*TEOl 


UNI«TEKX  MERPAR  ] 

HERPAR 


FIG.  A.  3  ROUTINES  OF  EXECUTIVE  PROGRAM  -  MAP 


Input  Data:  first  define  the  chart,  then  the  data 

PLEASE  INPUT  A  TITLE  (UP  TO  80  CHARS. ) 

> 

PLEASE  INPUT  YOUR  CHART  LIMITS 
IN  TERMS  OF  DEGREES  OF  LATITUDE  AND  LONGITUDE 
IN  THE  ORDER  TOP,  BOTTOM,  RIGHT,  LEFT 
(SOUTH  AND  WEST  NEGATIVE) 

> 

WHAT  IS  THE  CHART'S  STANDARD  PARALLEL  ? 

(DEGREES  AND  MINUTES  —  NORMALLY  THE  MID  LATITUDE) 

> 

AT  WHAT  NATURAL  SCALE? 

(F  WHERE  1/F  IS  THE  SCALE) 
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WHICH  SPHEROIDAL  CONSTANTS  DO  YOU  WISH  TO  USE? 
THOSE  FOR  THE  INTERNATIONAL  TYPE  1 
THOSE  FOR  THE  WGS  72  TYPE  2 

> 

DO  YOU  WANT  TO  PLOT  ANY  INTERMEDIATE 
PARALLELS  OF  LATITUDE? 

ANSWER  Y  FOR  YES  OR  N  FOR  NO 

> 


If  the  response  is  yes: 
HOW  MANY? 


PLEASE  LIST  THEM  AS  DEGREES  MINUTES 
(SOUTH  NEGATIVE) 

> 

TICK  MARKS  ON  THE  LATITUDE  SCALE? 

Y  FOR  YES  OR  N  FOR  NO 

> 


If  the  response  is  yes: 

TICK  INTERVAL  IN  DEGS.  AND  MINS? 

> 

DO  YOU  WANT  TO  PLOT  ANY  INTERMEDIATE 
MERIDIANS  OF  LONGITUDE? 

ANSWER  Y  FOR  YES  OR  N  FOR  NO 

> 


If  the  response  is  yes: 

HOW  MANY? 

> 

PLEASE  LIST  THEM  AS  DEGREES,  MINUTES 
(WEST  NEGATIVE) 

> 

TICK  MARKS  ON  THE  LONGITUDE  SCALE? 

Y  FOR  YES  OR  N  FOR  NO 

> 


If  the  response  is  yes  : 
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A. 3  Mean  Program 

Function:  To  compute  over  a  user-defined  temporal  and  spatial  field,  the 
mean,  maximum,  minimum,  standard  deviation,  gradient,  and  number  of 
observations  of  sound  speed  at  US  National  Oceanographic  Data  Center 
standard  depth  levels  <A.3>. 

Execution: 

3XQT  SMODS*SYSTEM. MEAN 
Subroutines:  (see  Fig.  A. 5) 


PIOBOR  ALMER  costa  herpar  lonlen  uni*tekx 


UNI*TEKX  BERPAR  |  [ 

BERPAR  [I  I  III  ||| 
HERPAR  CCP 


UNI*TEXX 


FIG.  A.  5  ROUTINES  OF  EXECUTIVE  PROGRAM  -  MEAN 


Input  Data: 

PLEASE  INPUT  THE  PLOT  WINDOW  AS  FOLLOWS 
MIN  AND  MAX  OF  THE  VELOCITY  SCALE  (M/S) 

MIN  AND  MAX  OF  THE  DEPTH  SCALE  (M) 

> 

PLEASE  INPUT  THE  CUT-OUT  DEPTH  AND  PLOT  SCALE  WHERE  THE  CUT-OUT  DEPTH 
IS  THE  MINIMUM  DEPTH  A  CAST  MUST  REACH  TO  BE  INCLUDED  IN  THE  ANALYSIS. 

THE  SCALE  WHERE  1.0  PRODUCES  A  PLOT  8  IN.  BY  5  IN. 

> 

The  CAMPO  routine  is  now  entered  to  define  the  data-field  limits  (see 
Sect.  2.5  of  main  text). 


Output  Data:  (see  Fig.  A. 6) 
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SOUND  VELOCITY  CHARACTERISTICS 
*********************** ************* 


DEPTH 

0 

NO 

14 

MIN 

1522.6 

MAX 

1536.6 

MEAN 

1529.7 

ST. DEV 
3.90 

GRAD 

10 

14 

1519.1 

1531.5 

1526.8 

3.54 

-28.43 

20 

14 

1515.0 

1524.9 

1520.1 

2.32 

-66.86 

30 

14 

1511.5 

1520.2 

1515.1 

2.46 

-50,43 

50 

14 

1506.3 

1513.7 

1509.1 

2.30 

-30.04 

75 

14 

1506.5 

1508.9 

1507,2 

.78 

-7.54 

100 

13 

1506.7 

1500.4 

1507.1 

.48 

-.37 

150 

12 

1507.5 

1508.2 

1507,9 

.25 

1.55 

200 

12 

150P.2 

1509.2 

1508,7 

.34 

1.55 

300 

12 

1509.9 

1511.1 

1510.3 

.36 

1.65 

400 

12 

1511.5 

1512.7 

1511.9 

.40 

1.61 

500 

11 

1513.1 

1514.2 

1513,5 

.37 

1.58 

FJ6.  A. 6 (a)  MEAN  -  PRINTOUT  OF  MEAN  PROFILE 


GRADIENT  M/S/ 1 00M 
□  V/DZ 

SO. SPEED  I  M/S  I  STANDARD  DEVIATION  I  M/S  ) 


FIG.  A. 6(b)  MEAN  -  PLOT  OF  MEAN  PROFILE,  STANDARD  DEVIATION  AND 

AVERAGE  GRADIENT 
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A. 4  CONTOUR  Program 

Function:  From  depth,  temperature,  salinity  or  sound  speed  data,  compute 
and  plot  contours  of  one  parameter  at  a  fixed  value  of  another  on  a 
Mercator  projection. 

Execution: 

qXQT  SMODS*SYSTEM. CONTOUR 
Subroutines:  (see  Fig.  A. 7) 


BATCH  HEAD  CHALIh  CAHPO  PARCHO  FILL  UNI*TEKX  NBGA  PAPER  CCP  COAST 


LONLEH  PLOBOR  ALMER  COSTA  HERPAR  UNl»TEKX 


UNI*TEKX 

FIG .  A.  7  ROUTINES  OF  EXECUTIVE  PROGRAM  -  CONTOUR 


Input  Data: 


The  CHALIM  routine  (see  Sect.  B.7  of  App.  B)  is  entered  first  to  define  the 
Meicator  chart,  followed  by: 


PLEASE  INPUT  THE  FIXED  AND  VARIABLE  PARAMETERS  TO  BE  CONTOURED 
BY  THE  FOLLOWING  KEY: 

1  FOR  DEPTH 

2  FOR  TEMPERA1URE 

3  FOR  SALINITY 

4  FOR  SOUND  SPEED 


> 

(e.g.  1  2  would  indicate  contours  of  temperature  at  a  fixed  depth  or  3  1 
would  indicate  contours  of  depth  of  a  fixed  value  of  salinity) 

INPUT  THE  LEVEL  AT  WHICH  TO  BE  CONTOURED 

(I.E.  OF  THE  FIXED  PARAMETER  E.G.  0  (ZERO)  METRES) 

> 
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At  this  point  the  CAMPO  subroutine  is  entered  to  define  the  data  field  (see 
Sect.  2.5  of  main  text)  followed  by  a  report  of  the  number  of  data  points 
which  have  been  extracted: 

DATA  EXTRACTION  COMPLETE 
(  )  DATA  POINTS  TO  BE  CONTOURED 

DO  YOU  WANT  THE  DATA  POINTS  TO  BE  PLOTTED? 

ANSWER  Y  FOR  YES 
OR  N  FOR  NO 

> 

PLEASE  INPUT 

THE  CONTOUR  INTERVAL 

CONTOUR  LABELLING  CODE  BY  : 

1  TO  3  (for  number  of  digits  required  after  decimal) 
or  -1  FOR  NO  DECIMAL  (I.E.  INTEGER  PART  ONLY) 
or  -3  FOR  NO  LABEL 
> 

THE  LINE-TYPE  CODE  BY  : 

1  THIN  LINE 

2  THICK  LINE 

3  DOTTED  LINE 

> 

At  this  point  the  plot  will  be  generated,  followed  by  a  request  for 
continuation  or  run  termination  as  follows  : 

PLEASE  INPUT  THE  DATA  CONTROL  KEYS  (1  OR  0) 

WHERE  KEY  1  **  NEW  TITLE 

KEY  2  **  NEW  DATA  FIELD  (TIME/SPACE) 

KEY  3  **  NEW  CONTOUR  LEVELS 
KEY  4  **  NEW  CHART  LIMITS 
KEY  5  **  NEW  CONTOUR  PARAMETERS 

> 

At  this  point  the  user  will  be  requested  to  input  only  those  data 
parameters  that  need  to  be  changed. 

A  response  of 

0  0  0  0  0 
will  terminate  the  run 


An  example  of  the  use  of  these  keys  is  as  follows: 

A  response  of 

10  10  1 

indicates  that  the  user  wishes  to  redefine  the  plot  title,  the  contour 
interval  and/or  line  type  and/or  label  type,  and  the  contour  parameters. 
In  such  a  case  only  those  questions  necessary  to  input  the  new  values  will 
be  asked.  The  remaining  data  fields  will  remain  the  same  as  defined  by  the 
previous  plot  generation. 


Output  Data:  (see  Fig.  A. 8) 
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8  CONTOUR  -  SEA  SURFACE  TEMPERATURE 
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A. 5  SECTION  Program 

Function  :  To  compute  and  plot  temperature,  salinity,  or  sound  speed  data 
as  a  function  of  depth  as  a  cross-section  between  two  input  geographical 
positions.  The  plot  displays  the  cross-section  in  terms  of  distance  in 
nautical  miles  from  the  initial  position  against  depth  in  metres,  either  as 
contour  lines  of  the  selected  parameter  or  as  data  points  or  both. 

Execution: 

SXQT  SMODS*SYSTEM. SECTION 
Subroutines:  (see  Fig.  A. 9) 


UNI*TEKX  PL030R  ALMER  COSTA  HERPAR  LONLEN 


UNI*TEKX  GENMAR 

NERPAR 

MERPAR 


CCP 


FIG.  A.  9  ROUTINES  OF  EXECUTIVE  PROGRAh  -  SECTION 


Input  Data: 

TITLE  OF  PLOT  (<  80  CHARS) 


At  this  point  the  CAMPO  subroutine  is  entered  to  define  the  data  field  (see 
Sect.  2.5  of  main  text) 

Followed  by: 

WHICH  PARAMETER  TO  SECTION? 

RESPOND  WITH  1  FOR  TEMPERATURE 

2  FOR  SALINITY 

3  FOR  SOUND  SPEED 

> 

TRACK  POSITIONS  IN  LAT  AND  LONG? 

(DEGS  AND  MINS  S  AND  W  NEGATIVE) 

> 

LANE  WIDTH  IN  N. MILES? (SWATH) 
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This  is  followed  by  plot  information: 

PLEASE  INPUT  MIN  AND  MAX  DEPTH  (METRES) 
RANGE  TICK  INTERVAL  (N.MI.) 
PLOT  SIZE  (X,Y,  IN  INCHES) 
CONTOUR  GRID  SIZE  (NX, NY) 

> 

PLOT  OF  DATA  POINTS  -  D 
CONTOURS  -  C 

OR  BOTH  -  B 

> 


If  response  is  D  or  B 

POINT  LABEL  CODE  AS  POSITIVE 
NUMBER  OF  FIGURES  AFTER  DECIMAL 
OR 

-  1  NO  DECIMAL  PART 

-  2  NO  LABEL 

-  3  NO  DATA  POINT 

> 


If  response  is  C  or  B 

CONTOUR  INTERVAL, LABEL  CODE, LINE  TYPE 

> 


At  this  point  the  plot  is  generated. 

Note:  The  data  field  selected  using  the  CAMPO  subroutine  defines  the 
catchment  area  of  the  data  base.  However,  this  executive  program  also 
checks  for  proximity  to  the  track  and  range  along  the  track.  This  range  is 
computed  as  shown  in  Fig.  A. 10,  by  projecting  the  position  of  the 
observation  on  to  the  track,  and  computing  the  distance  (d)  off  the  track 
and  the  range  (R)  from  the  initial  position. 

If  d  is  less  than  (SWATH/2)  the  observation  is  included.  After  all  the 
data  has  been  examined,  the  water  depths  of  included  data  are  sorted  by  the 
ORDINE  routine  (see  App.  B)  into  increasing  values  of  R  (range  along 
track),  such  that  the  bathymetric  profile  may  be  plotted.  The  parameter 
isovalues  are  extracted  and  the  contour  matrix  is  generated,  where  the 
abscissa  is  range  (R)  and  the  ordinate  is  depth. 


FIG.  A. 10  RANGE  CALCULATION  FOR  PROGRAM  SECTION 
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A. 6  SIGMAT  Program 

Function:  To  plot  or  printout,  within  user-defined  limits  of  temperature 
and  saMnity,  values  or  isopleths  of  sigma-t.  The  plot  may  be  superimposed 
with  user-selected  data  points  from  the  data  base. 

Execution: 

a  XQT  SMOD$*SYSTEM. SIGMAT 
Subroutines:  (see  Fig.  A. 12) 


UNI*TEKX  PLOBOR  AIMER  COSTA  MERPAR  LONLEN 


UNI*TEKX  MERPAR 

MERPAR 

I 

MERPAR 


CCP 

I 

UNI*TEKX 


FIG.  A.  12  ROUTINES  OF  EXECUTIVE  PROGRAM  -  SIGMAT 


Input  Data 

DO  YOU  WANT  OUTPUT  AS  A  LISTING  TYPE  L 
OR  PLOT  TYPE  P  ? 


If  response  is  L: 

PLEASE  INPUT  MIN  TEMP 
MAX  TEMP 
TEMP  INCREMENT 
MIN  SAL 
MAX  SAL 
SAL  INCREMENT 


> 
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if  response  is  P: 


PLEASE  INPUT  MIN  AND  MAX  TEMP.  VALUES 
MIN  AND  MAX  SAL.  VALUES 

> 

PLEASE  INPUT  LENGTH  (IN  INCHES) 

OF  TEMPERATURE  AXIS 
OF  SALINITY  AXIS 

> 

PLEASE  INPUT  SIGMAT  INTERVAL 

HEAVY  LINE  INTERVAL 
NO  OF  DIGITS  ON  CONTOUR  LABEL 
OR  -1  NO  DECIMAL  PART 
NO  LABEL 


DO  YOU  WISH  TO  SELECT  AND  PLOT  DATA  FROM  THE  O.C.  DATA  BASE 
Y  OR  N 


> 


If  the  response  is  Y  the  CAMPO  subroutine  is  entered  (see  Sect.  2.5  of 
the  main  text)  to  define  the  data  field. 

Output  Data:  (see  Fig.  A.  13) 


PLOT  OP  S1GMA-T 
SALINITY  l PP  r  ) 


30.0  40-0 


SALINITY  l  PPT  ) 


FIG.  A. 13 (a)  SIGMAT  -  ISOPLETHS  OF  SIGMA-t 


TEMPERATURE  ( C  ) 
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TEMPERATURE  t  C  ) 


**  ■5w^»;w:.!v’»‘r 


row .  ’■*/  /\  n 
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APPENDIX  B 
SUBROUTINE  LIBRARY 


B. 1  General 

Each  of  the  subroutines  used  by  the  system,  with  the  exception  of  those  of 
the  UNI*TEKX,  MFS  and  CCP  packages,  will  be  briefly  described  by  function, 
access  and,  where  applicable,  source  reference. 

They  may  also  be  used  independently  by  other  FIELDDATA  software,  and 
collected  by  the  EXEC  8  instruction 

LIB  SMODS*LIBRARY 

Data  are  passed  to  and  from  these  routines  either  by  formal  and  actual 
parameter  lists  or  by  access  to  one  or  more  of  a  series  of  named  COMMOM 
blocks. 

These  are: 

/C0NT/IHD(28) ,IDET(5,2500) ,NDET 

/C0NTR0/I0U , ITY.MODE ,MULT , 1ST , IPRNT , INDOUT , KBAT 

/FIELD/I0MCT(7) , IDO ,MSCT(25) , IMQ,IDQCT(25 , 100) , IDQ(25) , 

IYRCT(80) , IYR,MNCT(11) , IMN, ILAN ,NLAN(20) , IVES,NVES(20) ,MS 

/POSIZ/ALAT,ALO,AMN,ALN,NV,NY,TIM,XD 

/RAGGIO/KLAT , CLAM , KLON , CLOM , RAY , KAREA 

/DAREA/UP . DOWN , RIGHT , ALEFT 

/WINDOW/MIV ,MAV ,MID ,MAD , FAX , IPAR 

/BASE/SMP,W,FAC,XL,YL 

/CHART/ND , NM , SD , SM , ED , EM , WD , WM , SCD , SCM , NSC 
/GRID/IPALS,MERS, IPAL(20) ,PAM, (20) ,MER(20) ,AMER(20) , 

LAY , LAT I KD , LAT I  KM , LO Y , LOT I KD , LOT I KM 
/DP/IDP( 26) ,MT 

/P0INT/XP(2501) ,YP(2501j ,ZP(2501) ,MX,MY,FLV 
/BLANK/XB(1484) ,YB(1484) ,NB 

In  all  subroutine  descriptions,  standard  FORTRAN  naming  conventions  are 
applied:  the  first  letter  of  the  symbolic  name  of  a  formal  parameter  or 
constituent  parameters  of  a  common  block  is  any  of  the  six  letters 
I,J,K,L,M  or  N,  the  routine  requires  integer-type  data,  if  any  other 
character  is  used  then  real-type  data  are  required.  Automatic-type 
conversion  however  allows  for  conversion  of  integer  values  to  real  values 
without  affecting  a  computation,  whereas  the  converse  may  produce  a 
rounding-off  error  and  should  be  avoided. 

B.2  ALMER 

Computes  the  latitude  for  a  given  value  of  meridional  parts  on  the 
International  Speroid. 

CALL  ALMER  (A.B), 

whet  a  A  is  the  input  value  of  meridional  parts, 

B  is  the  computed  latitude. 

This  routine  requires  the  MERPAR  subroutine. 
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B.3  ALLSQ 

Computes  those  Marsden  Square  and  one  degree  square  references  in  and  on 
the  border  of  a  rectangle  defined  by  its  north,  south,  east  and  west 
geographical  limits 

CALL  ALLSQ 

This  routine  requires  input  data  passed  through  COMMON  blocks  FIELD  and 
DAREA  where  the  area  limits  are  input  in  DAREA  and  the  MSQ.DSQ  output  is 
put  into  FIELD. 

This  routine  requires  the  GENMAR  subroutine. 


B.4  BIGPIC 

Identifies  the  terminal  type  by  identification  of  the  terminal  site  ID. 
CALL  BIGPIC  (ID), 

where  ID  is  the  returned  terminal  type 

as  1  =  Batch  Run 

2  =  UNIVAC  DCT  500 

3  =  Tektronix  Model  4014 

4  =  Tektronix  Model  4010 

5  =  UNIVAC  UNISCOPE  100 

6  =  Tektronix  Model  4025. 

This  routine  requires  the  UNI*TEKX  subroutine  (not  described  here). 


B. 5  CALC 

Computes  sound  speed  <B.10>,  potential  temperature  <B.5>,  density  <B.7>, 
potential  density,  dynamic  height  anomaly  <B.7>,  specific  volume  anomaly, 
and  Brunt-Vaisala  frequency  <B.9>. 

CALL  CALC  (INST), 

where  INST  is  the  instrument  code  as  defined  in  Sect.  2.5  of  the  main  text. 
This  routine  requires: 

(a)  Common  blocks  CONT,  CONTRO,  POSIZ  and  FIELD 

(b)  routines  PRES,  WILSV,  SIG,  ROVO,  POT,  VBC  and  VBU 

(c)  UNI*TEKX. 

Output  is  to  the  unit  with  the  name  held  as  the  current  value  of  parameter 
IPRNT  of  block  CONTRO. 


B.b  CAMPO 

Oefines,  by  conversational  demand,  the  data  selection  field  limits  (see 
Sect.  2.5  of  the  main  text) 

CALL  CAMPO 
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This  routine  requires: 

(a)  COMMON  blocks  FIELD,  CHART,  GRID,  RAGGIO,  DAREA  and  BASE 

(b)  UNIMEKX. 

B. 7  CHALIM 

Defines,  by  conversational  demand,  the  limits,  meridians,  parallels, 
standard  parallel  and  natural  scale  of  a  Mercator  chart;  uses  these  data  to 
compute  the  chart-control  parameters  output  by  COMMON  block  BASE. 

CALL  CHALIM 
This  routine  requires: 

(a)  COMMON  blocks  CONTRO,  CHART,  GRID,  BASE,  FIELD 

(b)  routines  JANGLE,  LONLEN,  MERPAR 

(c)  UNI*TEKX. 

On  return,  BASE  will  hold: 

SMP  =  Meridional  parts  of  southern  border 
W  =  Longitude  in  minutes  of  western  border(West  negative) 

FAC  =  Mercator  Projection  scaling  factor  to  be  applied  to  all  N-S, 

E-W  differences 

XL,YL  =  the  chart  size  in  inches. 


B. 8  COAST 

Define,  by  conversational  demand,  the  input  parameters  for  routine 
COSTA.  CALL  CQAST 

This  routine  requires: 

(a)  COMMON  block  BASE 

(b)  routine  COSTA. 


C.9  CORD 

Converts  a  latitude  and  longitude  to  X/Y  coordinates  in  inches,  given  the 
necessary  control  data  (in  COMMON  block  BASE). 

CALL  CORD  (A,B,C,D) , 

where  A  =  Latitude  (degrees  ) 

B  =  Longitude  (degrees  ) 

C  =  X  value  (converted  longitude) 

D  =  Y  value  (converted  latitude). 

This  routine  requires: 

(a)  COMMON  block  BASE 

(b)  routine  MERPAR. 
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B. 10  COSTA 

Plots,  within  input  chart  limits,  the  coastline  selected  from  a  choice  of 
data  files  of  digitized  coordinates  (see  Sect.  3.2  of  the  main  text). 

CALL  COSTA  (I ,  J) , 

where  1  =  1  for  Mediterranean 

=  2  for  Southwestern  Approaches  to  the  English  Channel 
=  3  for  North  Atlantic. 

J  =  1  for  thin  line 
=  2  for  thick  1  ine 
=  3  for  dashed  line. 

This  routine  requires: 

(a)  COMMON  block  BASE 

(b)  routines  ERTRAN  (UNIVAC  FORTRAN  V  LIBRARY)  and  MERPAR 

(c)  UNI*TEKX. 


B. 11  DIST 

Compute  the  distance  and  course  between  two  geographic  positions  using  the 
principles  of  Mercator  and  Parallel  Sailing  <B.2>.  This  can  be  used  over 
only  relatively  short  distances,  (<  600  n.mi)  otherwise  the  great  circle 
distance  and  incremented  course  should  be  used. 

CALL  DIST  ( A, B,C,D, DIST .ANGLE)  , 

where  A,B  is  the  starting  position  in  degrees  of  latitude  and  longitude 

C,D  is  the  final  position  expressed  in  degrees  of  latitude  and 
longitude 

DIST  is  the  computed  distance  in  nautical  miles 

ANGLE  is  the  course  relative  to  north 
(i .  e.  +  or  -180°). 

This  routine  requires  the  MERPAR  subroutine. 


B.  12  EXPLAN 

Prints  out  an  explanation  (parameters  and  units)  of  the  computational 
output  of  the  CALC  subroutine. 

CALL  EXPLAN 

This  routine  requires  COMMON  block  CONTRO;  an  example  is  shown  in  Fig.  A. 2d 
of  App.  A. 

B. 13  FILL 

Uses  the  data  field  defined  by  CAMPO  and  the  data  provided  by  FINDBT,  to 
fill  the  data  matrix  contoured  by  the  CONTOUR  program  (see  App.  A. 4). 

CALL  FILL  (IP) 

where  IP  is  output  as  the  current  number  of  points  in  the  data  matrix. 


42 


SACLANTCEN  SM-151 


This  routine  requires 

(a)  COMMON  blocks  CONT,  BASE,  FIELD,  RAGGIO,  DAREA  and  POINT 

(b)  routines  FINOBT,  CORD 

(c)  MFS. 


B. 14  FINDBT 

The  interface  to  the  O.C.  Data  Base:  loads  into  core  a  profile  defined  by 
the  parameters  INST,  MSQ,  DSQ,  MON  and  NUMBER  <B.1>. 

CALL  FINDBT  (1,0, K,L,M,  (1),  (2)), 

where  I  =  Instrument  Code 

J  =  Marsden  Square  Number 
K  =  1°  Square  Number 
L  =  Month 

M  =  Consecutive  Element  No. , 
followed  by  two  error  returns: 

(1)  when  no  element  of  this  name  is  available  i.e.  no  such  data  exist 
within  this  I,J,K,L  combination 

(2)  when  no  file  of  this  name  is  available  i.e.  there  are  no  data 
within  the  combination  1,J,K,  (see  Sect.  3.3  of  Vol.  I  <B.1>). 

This  routine  requires: 

(a)  COMMON  block  CONT 

(b)  routine  STAND 

(c)  MFS. 

The  data,  if  retrieved  successfully,  are  loaded  into  COMMON  area  CONT. 


B.  15  GENMAR 

Computes  the  Marsden  Square  (MSQ)  and  1°  square  (DSQ)  reference  of  a 
geographical  position  (see  Figs. 2,  3  of  main  text). 

CALL  GENMAR  (A,B,K,L), 

where  A  the  latitude  in  degrees  \  .  . 

B  the  longitude  in  degrees  J  pu 

E  output 


B. 16  HEAD 

Prints  a  lint  iter-output  identification  for  any  batch  run. 
CALL  Hu..u 

This  routine  requires  COMMON  block  CONTRO. 
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B. 17  HEADER 

Prints  the  header  (or  dictionary)  data  of  the  data  base. 

CALL  HEADER  (INST, ICON), 

where  INST  is  the  Instrument  Code 

ICON  is  the  Consecutive  Number  of  the  profile  within  one 
INST/MSQ/DSQ/MONTH  group.  (Normally  a  loop  control  variable). 

Output  is  to  the  unit  with  the  current  value  of  parameters  IPRNT. 

This  routine  requires  COMMON  blocks  CONT,  CONTRO,  FIELD  and  POSIZ. 

B. 18  INCHEK 

Checks  whether  a  profile  is  within  certain  of  the  data  fields  selected  by 
S/R  CAMPO. 

CALL  INCHEK  ($(1)), 

where  control  is  returned  normally  if  a  profile  lies  within  the  fields  of 
spatial  definition,  year,  country,  and  ship;  in  addition  the  routine  checks 
the  "  doubtful  "  flag  for  inclusion  or  otherwise.  If  any  one  of  these 
field  limits  is  violated,  control  is  returned  through  label  $(1). 

This  routine  requires  COMMON  blocks  CONT,  CONTRO,  FIELD,  RAGGIO  and  DAREA. 


B.  19  JANGLE 

Rings  the  terminal  bell. 

CALL  JANGLE  (N), 

where  N  =  number  of  rings  required. 
This  routine  requires  UNI*TEKX. 


B. 20  LIMITS 

Allows  the  user  to  delineate  a  sub-area  of  the  Atlantic  Ocean  by  inter¬ 
active  use  of  the  Tektronix  cursor.  Returns  the  limits  of  the  selected 
area,  together  with  a  computed  natural  scale  and  mid- latitude  required  to 
plot  the  area  within  the  limits  of  the  Tektronix  small  screen  (i.e.  models 
4002,  4010,  4012). 

CALL  LIMITS 
This  routine  requires 

(a)  COMMON  blocks  BASE,  CHART,  GRID,  DAREA 

(b)  routines  LONLEN,  MERPAR,  PL0B0R,  COSTA,  ALMER 

(c)  UNIMEKX. 

On  entry  into  this  routine  the  following  message  is  output: 

A  MAP  OF  THE  ATLANTIC  OCEAN  WILL  NOW  BE  PLOTTED  FOR  YOU 
TO  INTERACTIVELY. 

INPUT  YOUR  DATA  AREA. 

WHEN  THE  CURSOR  IS  ENABLED,  INPUT  THE  SOUTH-WEST  AND  NORTH-EAST 
CORNERS. 
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At  this  point,  the  coastline  is  plotted  and  the  cursor  enabled.  If  the 
input  positions  for  the  corners  are  not  of  the  right  order  of  magnitude  or 
are  input  in  an  incorrect  order,  an  error  message  is  output  and  the  input 
is  re~sol icited. 


B. 21  LIST 


Prints  a  station  of  the  data  base,  with  the  option  to  print  observed  and/or 
standard  depth  levels.  Does  not  produce  a  character  plot  (see  S/R  LPLOT). 

CALL  LIST  (INST,  ICON), 

where  INST  =  Instrument  Code 

ICON  =  Consecutive  number  of  the  profile  within  one 
INST /MSQ/DSQ/MONTH  group 

This  routine  requires 

(a)  COMMON  blocks  CONT,  CONTRO.  FIELD  and  POSIZ 

(b)  routine  TESTA. 

Output  is  to  the  unit  assigned  with  the  current  value  of  variable  IPRNT  of 
/CONTRO/. 


0. 22  LONLEN 

Computes  the  length  in  metres  of  1'  longitude  at  a  given  latitude; 
including  an  option  of  reference  ellipsoid  (see  Sect.  3.3  of  main  text). 

CALL  LONLEN  (A,Z,I), 

where  A  =  latitude  in  degrees 

Z  =  computed  length  of  1'  longitude  (metres)  at  lat  A. 

1  =  1  for  International  Spheroid 
2  for  WGS72. 


B.  23  LPLOT 

Prints  a  station  of  the  data  base  (as  S/R  LIST),  including  a  character  plot 
of  each  component  parameter  profile  (Fig.  A. lb  of  App.  A) 

CALL  LPLOT  (INST, ICON) , 

where  INST  =  Instrument  Code  (I) 

ICON  =  Consecutive  number  of  the  profile  within  one 
INST/MSQ/DSQ/MONTH  group 


B. 24  MERPAR 

Computes  the  meridional  parts  of  a  given  latitude,  with  the  option  of 
reference  ellipsoid. 

CALL  MERPAR  (A.Z.I), 
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where  A  =  Latitude  in  degrees 
Z  =  Meridional  parts 
1  =  1=  International  spheroid 
2  =  WGS72. 


B.  25  MUST 

Prints  the  results  of  the  executive  MEAN  program  (see  App.  A. 3). 

CALL  MLIST  (B,IA) 

where  B  =  output  array  dimensioned  for  five  output  parameters  at  up  to  26 
standard  levels. 

IA  =  Array  of  standard  depth  levels. 

This  routine  requires  COMMON  blocks  CONTRO,  DP. 

Output  is  to  the  unit  assigned  with  the  current  value  of  parameter  IPRNT  of 
/CONTRO/;  an  example  is  shown  as  Fig.  A. 6  of  App.  A. 


B. 26  MPLOTS 

Plots  the  output  of  program  MEAN  (see  App.  A. 3). 

CALL  MPLOTS  (TITLE) 

where  TITLE  is  a  character  string  shorter  than  81  FIELDDATA  characters. 

This  routine  requires  common  blocks  WINDOW  and  DP. 

An  example  of  the  output  is  shown  as  Fig.  A. 56  of  App.  A. 


B. 27  NBGA 

Checks  whether  a  variable  has  the  value  of  either  Y  or  N  (for  Yes  or  No)  or 
whether  it  has  another  value. 

CALL  NBGA  (CHAR,  $) 

Returns  normally  if  CHAR  has  the  value  Y  or  N;  returns  through  label  $  if 
any  other  value.  It  is  designed  essentially  to  "  idiot-proof  "  the 
conversational  software. 


B. 28  NBGN 

Checks  whether  a  numerical  value  lies  within  a  certain  range. 

CALL  NBGN  (ALOW, HIGH, VAL,  $) 

where  ALOW  is  lowest  permissible  value 
HIGH  is  highest  permissible  value 
VAL  is  the  value  being  tested. 

The  routine  returns  control  normally  if: 

ALOW  <  VAL  <  HIGH 

otherwise  it  returns  through  label  $. 
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B. 29  NDIG 

Determines  the  number  of  digits  in  an  integer  value;  used  in  character 
spacing  on  graphic  displays.  This  is  a  FORTRAN  Function  and  is  therefore 
accessed  by: 

(Variable)  =  NDIG(I), 

where  I  =  the  value  being  plotted. 

B. 30  ORPINE 

Sorts  a  variable  number  of  paired  integer  and  real  values  into  ascending 
order  of  magnitude  of  the  real  values,  using  the  "bubble  "  sorting  method 
(Knuth  <B.3>). 

CALL  ORPINE  (A.I.N), 

where  A  is  a  real  array 

I  is  an  integer  array 
N  is  the  number  of  pairs  of  values 

(i.e.  arrays  are  dimensional  A(N),I(N)). 

B.  31  OUTPUT 

Sets  up  the  control  variables  to  direct  the  output  of  executive  program 
DISPLAY.  The  conversational  demands  are  those  explained  in  App.  A.l. 

CALL  OUTPUT 

This  routine  requires 

(a)  COMMON  blocks  C0NTR0,  WINDOW,  FIELD 

(b)  UNIHEKX. 


B. 32  PAPER 

A  Calcomp  960  drum  plotter  is  employed  at  SACLANTCEN.  A  real-time  program 
schedules  plots  generated  by  user  programs  on  to  the  plotter  in  a  range  of 

fixed  paper  sizes  <B.4>.  Given  a  plot  size  in  inches,  this  routine 

computes  the  smallest  permissible  sheet  size  and  its  orientation,  onto 
which  the  plot  will  fit. 

CALL  PAPER  (X,Y,IP,  $), 

where  : 

X,Y  are  the  abscissa  and  ordinate  lengths  in  inches 
IP  is  the  plot  size  in  the  range  1  to  12 

$  is  an  error  return  label  if  the  plot  size  is  greater  than  the 

maximum  possible  sheet  size. 


B. 33  PARCHO 

Solicits  information  of  required  contour  parameter  and  level  for  the 
CONTOUR  executive  program. 

CALL  PARCHO 

This  routine  requires  COMMON  block  POINT. 
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B. 34  PLOBOR 

Plots  a  Mercator  chart  border,  with  intermediate  lines  of  latitude  and 
longitude  labelled  in  terms  of  degrees  and  minutes  E  or  W,  N  or  S.  In 
addition  plots  tick  marks  on  either  or  both  scales  as  required  (see 
Fig.  A. 4a  of  App.  A) 

CALL  PLOBOR 

This  routine  requires: 

(a)  COMMON  blocks  CHART,  BASE,  GRID 

(b)  routines  NDIG,  MERPAR 

(c)  UNI*TEKX. 


B. 35  POT 

Computes  potential  temperature  from  in-situ  pressure,  temperature  and 
salinity,  using  Fofonoff's  equation  <B.5>. 

V  =  POT  (P,T,$), 

where  P  =  pressure  in  decibars 
T  =  temperature  in  °C 
S  =  sal inity  %0 

and  V  =  potential  temperature  °C. 


B. 36  PRES 


Converts  depth  to  pressure  at  a  latitude,  using  Leroy's  1968  <B.6> 
simplified  formula. 

V  =  PRES  (D,  ALAT), 
where  D  =  depth  in  metres  (F) 

ALAT  =  latitude  in  degrees 
and  V  =  pressure  in  kg/cm2  (F). 


B. 37  QUAD 

Given  two  positions,  computes  the  course  from  the  first  to  the  second  as 

(a)  a  signed  angle  from  0°  to  90°  whi  h  by  inspection  may  be 
converted  to  cardinal  notation,  (see  Fig  B.l) 

(b)  the  course  expressed  in  three-figure  notation,  e.g.  if  a  =45° 

the  cardinal  notation  is  S65°W  and  the  three-figure  notation  is 
245°;  conversely  if  a  =  -65°  the  cardinal  notation  is  S65°E  and 
the  three-figure-notation  is  115°. 

CALL  QUAD  (Al,  01  ,A2,  02,  C90,  C360), 

where  Al,01  are  the  latitude  and  longitude  of  position  1 
A2,02  are  the  latitude  and  longitude  of  position  2 
(S  and  W  negative) 

C90  is  the  signed  cardinal  angle  from  PI  to  P2 
C360  is  the  360°  angle. 
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FIG.  B.l  CARDINAL  ANGLE  TO  THREE-FIGURE  NOTATION 
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B. 38  RADIAL 

Given  a  circle  of  radius  R  on  the  earth's  surface  centred  at  position  ATL, 
OTL,  computes  the  geographical  limits  of  the  enclosing  square  (see  Fig.  5 
of  main  text). 

CALL  RADIAL  (ATL,  OTL,  RAD), 

where  ATL,  OTL  is  the  position  of  the  centre  (latitude,  longitude) 

RAD  is  the  radius  (n.mi). 

Computed  results  are  stored  in  COMMON  area  DAREA,  which  is  therefore 
required. 


B. 39  ROVO 

Computes  density  and  specific  volume  from  sigma-t, pressure,  temperature, 
and  salinity  <B.7>. 

CALL  ROVO  (SG,  P,  T,  S,  RHO,  SVO), 

where  SG  is  sigma-t  in  (g/cm3  -  1)/1000 
P  is  pressure  in  decibars 
T  is  temperature  in  °C 
S  salinity  in  %0 
RHO  is  density  in  g/cm3 
SVO  is  specific  volume  (1/RHO). 

B. 40  SIG 

Computes  sigma-t  from  temperature  and  salinity  <B.7>. 

V  =  SIG  (T,  S), 

where  T  =  temperature  in  °C 
S  =  salinity  in  % 

giving  V  the  sigma-t  in  (g/cm3-l)/1000. 


B.44  STAND 

Interpolates  significant  point  data  to  US  National  Oceanographic  Data 
Center  standard  depths  <B.8>  using  linear  interpolation. 

In  addition,  each  observed  depth  is  labelled  with  a  3,  each  interpolated 
depth  is  labelled  with  a  6. 

CALL  STAND  (INS), 

where  INS  is  the  Instrument  Type  Code  (Sect.  2.5  of  the  main  text). 

This  routine  requires  COMMON  block  CONT. 

The  computed  results  are  merged  with  the  input  observed  profile  in  area 
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B.42  TESTA 

Re-formats  the  data-base  dictionary  into  correct  units. 

CALL  TESTA(INST) , 

where  INST  is  the  Instrument  Type  Code  (Sect.  2.5  of  main  text). 

This  routine  requires  COMMON  blocks  CONT  and  POSIZ.  The  reformatted 
variables  are  output  in  area  POSIZ. 


B. 43  1PL0T 

Plots  single  or  multiple  profiles  of  temperature,  salinity,  or  sound  speed 
as  a  function  of  depth. 

CALL  TP LOT  (INST) 

where  INST  is  the  Instrument  Type  Code  (Sect.  2.5  of  main  text). 

This  routine  requires 

(a)  COMMON  blocks  C0NT.C0NTR0  and  WINDOW 

(b)  UNI*TEKX. 

The  plot  is  self scaled  within  the  input  values  held  in  WINDOW,  and  the  axes 
labelled  accordingly.  Single  profile  plotting  initiates  and  closes  a  plot 
file  (see  Sect.  4.2  of  main  text)  each  time  the  routine  is  called,  whereas 
multiple  profile  plotting  only  initiates  on  the  first  entry  and  closes  on 
the  last.  Examples  of  the  output  are  shown  as  Figs.  A.2e  and  A.2f  of 
App.  A. 


B. 44  UKLIST 

A  modified  version  of  routine  LIST  to  output  United  Kingdom  Hydrographic 
Office  formatted  meteorological  data  <see  Sect.  3.3  of  Vol.  I  <B.1>). 

CALL  UKLIST  (INST, ICON), 

where  INST  =  Instrument  Type  Code  (Sect.  2.5  of  the  main  text) 

ICON  =  Consecutive  number  of  the  profile  within  one 

INST/MSQ/DSQ/MONTH  group  (normally  a  loop  control  variable). 

This  routine  requires  : 

(a)  COMMON  blocks  CONTRO,  CONT  and  POSIZ 

(b)  routine  TESTA. 

Output  is  written  to  the  data  file  or  unit  having  the  name  assigned  with 
the  current  value  of  variable  IPRNT  of  CONTRO;  an  example  is  shown  as  Fig. 
A. 2c  of  App.  A. 
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B. 45  VBC 

A  function  to  compute  Brunt-Vaisala  frequency  <B.9>. 

V  =  VBU  (DG,G) 

where  DG  =  density  gradient 

G  =  gravity, 

giving  V  in  cycles/hour. 

B. 46  WILSV 

Computes  the  speed  of  sound  in  water  using  Wilson's  2nd  equation  <B.1C>. 

V  =  WILSV  (T,S,P), 

where  T  =  temperature  (°C) 

S  =  salinity  (%) 

P  =  pressure  kg/cm2 , 

giving  V  =  sound  speed  (m/s). 
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COMNAVSOUTH  1 

COMSTRIKFORSOUTH  1 

COMEDCENT  1 

COMMARAIRMED  1 

CINCHAN  1 


Copies 

SCNR  FOR  SACLANTCEN 

SCNR  Belgium  1 

SCNR  Canada  1 

SCNR  Denmark  1 

SCNR  Germany  1 

SCNR  Greece  1 

SCNR  Italy  1 


SCNR  Netherlands  1 
SCNR  Norway  1 
SCNR  Portugal  1 
SCNR  Turkey  1 
SCNR  U.K.  1 
SCNR  U.S.  2 
SECGEN  Rep.  SCNR  1 
NAMILCOM  Rep.  SCNR  1 


NATIONAL  LIAISON  OFFICERS 


NLO  Canada  1 
NLO  Denmark  1 
NLO  Germany  1 
NLO  Italy  1 
NLO  U.K.  1 
NLO  U.S.  1 

NLR  TO  SACLANT 

NLR  Belgium  1 
NLR  Canada  1 
NLR  Denmark  1 
NLR  Germany  1 
NLR  Greece  1 
NLR  Italy  1 
NLR  Netherlands  1 
NLR  Norway  1 
NLR  Portugal  1 
NLR  Turkey  1 
NLR  UK  1 
NLR  US  1 

Total  initial  distribution  236 
SACLANTCEN  Library  10 
Stock 

Total  number  of  copies  280 


